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We show that Larkin-Ovchinnikov (LO) states with modulated supernuid order parameters have 
a considerably larger range of stability in a lattice than in the continuum. We obtain the phase 
diagram for the 3D cubic attractive Hubbard model with an unequal population of up and down 
fermions using the Bogoliubov-de Gennes fully self-consistent method. We find a strong modulation 
of the local polarization that should provide a distinct signature for detection of the LO phase. 
The shell structure in the presence of a trap generates singularities in the density at the phase 
boundaries and provide additional evidence for the LO phase. Depending on specific parameters, 
the LO ground state occurs over a large range of population imbalance, involving 80% of the atoms 
in the trap, and can exist up to an entropy s ~ 0.5fcs per fermion. 



An imbalanced population of fermions with two hy- 
pcrfine states and interacting via attractive interactions 
offers the exciting possibility of observing superfluidity 
with a spatially modulated order parameter. For small 
imbalance the ground state is a BCS/BEC supernuid 
state with paired fermions, but for large imbalance the 
ground state is a polarized Fermi liquid, [l], Q At inter- 
mediate polarizations, mean- field theories predict Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) states with a spa- 
tially modulated supernuid order parameter 0, 0| that is 
a compromise between pairing and polarization. There 
is controversy over whether the FFLO state exists. ID 
systems only allow a quasi-long-range-ordered version of 
FFLO, whereas in 2D and 3D continua, FFLO only oc- 
cupies a sliver of the phase diagram and is vulnerable to 
fluctuations. So far, ordered FFLO has not been ob- 
served except in some reports on layered organic and 
heavy-fermion superconductors [Bj]. Both BCS and po- 
larized states have been observed in imbalanced cold 
fermionic gases 0, 0], but the LO phase has so far re- 
mained elusive. 

In this paper we study the full phase diagram of the 
cubic lattice Hubbard model. We use approaches based 
on variational mean-field theory in six channels, which 
includes Bogoliubov-de Gennes (BdG) and Hartree cor- 
rections. Our main results are: (1) The LO phase oc- 
cupies a large region of the phase diagram between the 
BCS and polarized Fermi liquid (2FL) phases. (2) With 
increasing field (or imbalance), the BCS state becomes 
unstable to an LO phase consisting of domain walls at 
which the order parameter changes sign. The polariza- 
tion is confined to these domain walls. At higher fields 
the domain wall structure evolves into a sinusoidal vari- 
ation of the order parameter accompanied by a polariza- 
tion variation at twice the wavevector. We suggest that 
the most promising way to detect the LO phase is to fo- 
cus on this spatial variation of the polarization. (3) The 
momentum distribution functions n a (k) in the LO phase 
shows distinct features that break the lattice symmetry, 
such as Fermi arcs, Fermi pockets, and blocking regions, 
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FIG. 1: Schematic depiction of an imbalanced fermion gas 
in an optical lattice in a cigar-shaped trap, corresponding to 
slice (i) in Fig. 2(b) There are regions of strong LO order con- 
sisting of domain walls (DWs), in which the excess majority- 
spins fermions are confined. The polarization in each DW can 
be about 0.25/a 2 fermions per unit area, the DW thickness 
about 3a, and the DW spacing about 6a (where a is the lattice 
constant). 



unlike in the BCS/BEC state. (4) Depending on parame- 
ters, the LO phase can exist below an entropy s ~ O.bks- 
(5) In an optical lattice in a shallow trap with appropri- 
ate parameters, LDA predicts that most (> 80%) of the 
atoms participate in the LO phase! See Fig. [1] 

Model and methods: The Hubbard Hamiltonian is 
given by 
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where £k = —2{cosk x + cosk y + cosfc 2 ) is the disper- 



sion relation on the cubic lattice for nearest-neigbor hop- 
ping, a = ±1 labels (hyperfhie) spin states, n Ta = c\ a c T(J 
are number operators, \x„ = fj, + ah are the chemical 
potentials for the two spin species, and U is the local 
pairwise Hubbard interaction. The hopping t is the unit 
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TABLE I: Effect of the "Lieb-Mattis transformation" (LMT) 
on the parameters, observables, and phases of the Hubbard 
model on a bipartite lattice. The LMT is a particle-hole 
transformation on the down spins with a it phase shift for 
the B sublattice: ca{ — > c A ,, cb{ — > — c B ,. It relates mag- 
netic phases of the repulsive model to paired/density-ordered 
phases of the attractive model. A Z-7T-SDW has a spin den- 
sity modulation in the z direction at the antiferromagnetic 
wavevector (-7T, 7r, tt), whereas an i-q-SDW has moments in 
the x direction at some general wavevector q; the latter maps 
to a LO (Larkin-Ovchinnikov) state, as illustrated in Fig. [3] 



of energy. We use the convention that repulsive U is 
positive. We find it convenient to work in terms of 
the average chemical potential p and the Zeeman field 
h = 2(/Xf — hi). The observables of interest are the den- 
sity n r = ^(n r i + n r i), imbalance m r = \{ n r] — n ri): 
and pairing density F r = (c r jC r |) = ^j. Since H has 
been written in terms of symmetrized densities n a — \ 
measured with respect to half-filling, the thermodynamic 
quantities and phase diagram are invariant under sign 
changes of p and h. Furthermore, the "Lieb-Mattis" 
transformation (LMT) relates the repulsive and Hubbard 
model as described in Table Q] 

Our calculations are based on Tr p In p variational 
mean-field theory. The Hubbard U is approximated by 
67V potentials, where N is the number of sites: the local 
chemical potentials /i" 1 ', Zeeman fields h* ot , and com- 
plex pairing potentials A r . This is equivalent to full BdG 
with Hartrce corrections. The trial Hamiltonian involves 
a 47V x AN matrix H: 



FIG. 2: Phase diagrams of the cubic lattice Hubbard model as 
functions of average chemical potential p and Zeeman field h. 
At large negative (positive) p the system is an insulator with 
every site empty (doubly occupied). At large h it is a fully- 
polarized insulator. 1FL represents a fully polarized Fermi 
liquid (half-metal). 2FL represents a two-component Fermi 
liquid. BCS and LO represent superfluid states. The black 
curves indicate Lifshitz transition boundaries corresponding 
to changes in Fermi surface topology. The red curves are 
slivers of a CDW phase that only occur at extremely low 
temperature, and are irrelevant to experiments. The dashed 
blue lines represent slices through the phase diagram corre- 
sponding to two different combinations of overall polarization 
and trapping potential; each blue dot indicates the chemical 
potential at the center of the trap. Slice (i) corresponds to 
a large LO core surrounded by 2FL and 1FL shells, as illus- 
trated in Fig. [T] Slice (ii) corresponds to an insulating core 
with two fermions per site, surrounded by seven shells. 



trix p and thus to a variational free energy "Tr p In p" 
^var = -f tr In (2 cosh + [ l 
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where the "pairing density' F T = (c r |C r |), charge density 
n T , and spin density m r are elements of the matrix Green 
function G = — i tanh ^ . Differentiating f2 var with re- 
spect to the variational parameters leads to the BdG and 
Hartree self-consistency conditions, and computing f2 var 
itself is necessary to distinguish between various possible 
ground states. We use simplified methods (e.g., instabil- 
ity analysis) where applicable. 

Phase diagrams: The ground state phase diagram has 
three parameters: p, h, and U. It is symmetrical un- 
der sign changes of h and p. Fig. 2(a) shows the non- 
interacting (U = 0) phase diagram. At finite attractive 
\U\, two superfluid phases (BCS and LO) appear in the 
phase diagram. The size of the LO region grows with \U\ 
until |£/| « 5t. At large \U\, mean-field theory is affected 
by fluctuations. Therefore we will focus on U = —6t 
(see Fig. |2(b)[ ). Note that U u = -7.91355i is the cou- 
This trial Hamiltonian corresponds to a trial density ma- pling where a bound state first appears in the lattice 
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(the analogue of unitarity in the continuum), while the 
bandwidth is I2t. Hartree corrections cause the phase 
diagram to be squashed in the /^.-direction and elongated 
in the /i-direction. The main encouraging observation is 
that between the BCS state at h = and the polarized 
FL state at moderate h, there is a sizeable region where 
the ground state is an LO state. 
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(c)Weak texture 
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FIG. 3: (a) Weak LO state with sinusoidal pairing density 
F r ~ cos q • r accompanied by CDW and SDW at wavevector 
2q. (b) Strong LO state just above the critical field for domain 
wall penetration. At each domain wall, the order parameter 
changes sign, and the polarization is finite due to occupation 
of Andreev bound states, (c) Weakly textured state, in which 
the z-components of moments are nearly uniform and the x 
components oscillate at wavevector (it, n, n) — q. (d) Strongly 
textured state (not to scale). The red circles represent the 
hole density with respect to half filling, which is concentrated 
at the domain walls of the antiferromagnet. Panels (a) and 
(c) are related to (b) and (d) by the LMT. 

Nature of the LO state: Figure [3] shows examples of 
weak and strong LO ground states in real space, and the 
textures in the repulsive Hubbard model that they corre- 
spond to under the LMT. At high fields just below h^o, 
there is a weak LO state characterized by a sinusoidal 
order parameter with a q- vector related to the difference 
between the Fermi wave vectors of the majority and mi- 
nority components. At lower fields this crosses over to 
a strong LO state consisting of BCS-like regions sepa- 
rated by domain wallsjH, [Iq|. The LO region actually 
includes a rich variety of patterns with different q (verti- 
cal stripes, diagonal stripes, 2D modulations) depending 
on commcnsurability effects. 

Our variational calculations find that FFLO ground 
states are always LO states that have a real order pa- 
rameter and break translational symmetry. These states 
have a pairing energy that can be 50 times larger than 
the FF states studied in Ref. @ that have a complex or- 
der parameter A ~ e lq r and break time-reversal symme- 
try. This observation is consistent with previous results 
for the continuum. [il, fl~ol] According to the LMT, an LO 
state maps to a coplanar spin texture (spins in xz-planc) , 
whereas an FF state maps to a non-coplanar "helical" 
texture. Thus, for the repulsive Hubbard model, copla- 
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FIG. 4: Momentum distributions of up- and down- fermions 
in an LO state with wavevector ±q = (±2^,0,0). rcj.(k) 
shows the effect of "blocking " 0] by the sea of up spins 
n-f (k) shifted by ±q, and vice versa. The steps in the orange 
curve correspond to remnants of the up-spin Fermi surface. 
The pattern breaks fourfold symmetry but preserves inver- 
sion symmetry, in contrast to Fig. 3 of Ref. @, where the e* q r 
ansatz breaks inversion symmetry. 



nar textures are more favorable. 

The pairing of up and down fermions belonging to un- 
equal sized Fermi surfaces leads to complicated features 
in the momentum distribution function n a (k), such as 
Fermi arcs, Fermi pockets, and blocking regions. Many 
patterns are possible depending on the parameters of the 
LO state. Figure 2] shows an example of the momen- 
tum distributions, in which selective pairing produces 
non-monotonic behavior. The most robust feature is the 
breaking of the lattice symmetry. In experiments this 
effect may be complicated by twinning due to trap ge- 
ometry. 

Shell structure in a trap: We now apply the above 
results for n^r) and n a {k) to optical lattices in 
traps, within the local density approximation (LDA), 
which should be applicable to shallow traps with many 
fermions. In LDA, the local phase is assumed to be de- 
termined by the local chemical potential, /i(r) = fiQ — 
VtrapM- This predicts shell structures corresponding to 



horizontal slices through the phase diagram (Fig. 2(b)); 
one such structure is depicted in Fig. [TJ Spherical traps 
may cause twinning between LO states of different orien- 
tations, whereas a cigar-shaped trap helps align domain 
walls perpendicular to the long axis. In principle, 2FL- 
LO (and/or LO-BCS) transitions show up as kinks in the 
density profiles n^(r) and n-i(j) with appropriate critical 
exponents; however, whether these kinks are observable 
depends on parameters and experimental resolution. 

Entropy for observing LO states: We predict that LO 
phases should be possible to observe at temperatures (or 
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(a)Phasc diagram for U = -6t, h = 1.5* (b)fi = 0.5* (BCS g.s.) (c)/i = 0.8* (LO g.s.) 



FIG. 5: (a) A finite-temperature phase diagram, (b) and (c): Entropy s (per site), average polarization m, and pairing 
F for /i = 0.5* and [i = 0.8* respectively, corresponding to dashed lines in (a). The BCS phase has m = at T = 0; m 
becomes non-zero at finite T due to polarized quasiparticle excitations. In contrast, the LO phase has m > even at T — 0. 
Note the different horizontal and vertical scales. Crosses indicate critical temperatures and entropies (sbcs ~ 0.69fcs and 
slo = 0.29fcs); s — ks In 4 = 1.386&B is the entropy of the Hubbard model at T = oo. 



entropies) that arc accessible to experiments. Within 
MFT we have found LO states up to Tlo = 0.6*. (slo — 
0.5fcs), for the parameters U = — 6t, fi = 0.254, h = 2.25*. 
This is not much lower than the critical temperature 
for the BCS phase, Tbcs ~ 1-1* (sbcs ~ O.Sfcs), at 
U = -6t, fjt = 0.25*, h = 0. Ref. [HI implies that the 
critical entropy of the cubic lattice BCS phase in the 
\U\ 12t regime is only reduced by about a factor of 
2 from its mean-field value (fcgln2). Since our results 
are at medium coupling, \U\ = 6t, we expect our MFT 
estimate of sbcs to be even closer to the true value. Fig- 
ure [S] compares the effect of temperature on BCS and LO 
phases. 

In conclusion, we find that for the cubic lattice within 
fully self-consistent mean-field theory, LO states occur 
over an enhanced region of the phase diagram, as com- 
pared to the continuum. This suggests that unbal- 
anced ultracold fermion systems in optical lattices should 
readily exhibit LO ground states, which could be de- 
tectable by virtue of the accompanying polarization os- 
cillations. Based on our calculations, we find that for 
N ~ 10 5 fermions with an overall polarization P ~ 0.37 
at coupling U — ~6t, about 83% of the atoms are in 
the LO phase. The polarization in each domain wall 



P 



DW 



E 



7t T (r)-n i (r) 



for a strong LO state such 
as in Fig. |3(b)| is about 30%; the polarization between 
domain walls is practically zero, giving a large contrast. 
The spacing between domain walls can be of order 10a, 
where a is the optical lattice constant. In optical lattice 
experiments with a ~ lfxxn, which implies a domain wall 
spacing of about 10/im. 

We explicitly describe the connection between LO 
phases in the attractive Hubbard model and spin textures 
in the repulsive Hubbard model via the LMT (a mapping 
that does not exist in the continuum). Thus, by chang- 
ing the sign of U, experiments in traps can effectively 
measure slices through the phase diagram in Fig. |2(b)| 



in both horizontal and vertical directions. We also point 
out that the strongly attractive Hubbard model at half- 
filling and weak imbalance will have a tendency towards 
exotic d-wave magnetism described by order parameters 

such as ^cj^Cj^y ~ cosfc^ — cosfcj,, as a consequence of 
applying the LMT to the repulsive Hubbard model with 
its d-wavc pairing tendencies. 

Given that the enhanced LO regions in the cubic lat- 
tice occurs because of nesting, we expect that anisotropy 
will improve nesting and further enhance the LO re- 
gion. Other authors have studied 2D arrays of ID tubes 
3 HI; ^ remains to be shown whether anisotropic lat- 
tices or coupled tubes are more favorable for LO. It will 
be important to include quantum and thermal phase fluc- 
tuations in reduced dimensions to get accurate estimates 
of phase boundaries. 
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